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GN ' Abstract 

Qh| We present a simple proof of the absence of Bose-Einstein condensation of 

P-i! a relativistic boson gas, in any finite local magnetic field in less than five 

dimensions. We show that the relativistic charged boson gas exhibit a genuine 
Meissner-Ochsenfeld effect of the Schafroth form at fixed supercritical density. 
As in the well-known non-relativistic case, this total expulsion of a magnetic 
field is caused by the condensation of the Bose gas at vanishing magnetic field. 
The result is discussed in the context of kaon condensation in neutron stars. 



1 Email address: elmfors@nordita.dk. 
2 Email address: tfepl@fy.chalmers.se. 



3 Email address: tfedp@fy.chalmers.se. 

4 Email address: tfebss@fy.chalmers.se. Research supported by the Swedish National Research Council 
under contract no. 8244-311. 



In this letter we study some aspects of an ideal charged Bose gas at finite temperature T 
and chemical potential \x in presence of a static uniform magnetic field. Magnetic fields 
B associated with compact astrophysical objects may range between B = O(10 4 ) T for 
magnetic white dwarfs to B = O(10 10 ) T for supernovae PJ. As a reference we recall that 
the characteristic magnetic field in QED is m 2 e /e = O(10 9 ) T. For such systems one may 
expect that also thermal effects are of importance. In the absence of a magnetic field it 
is known that the relativistic charged boson gas exhibits a Bose-Einstein condensation 
(see e.g. Ref. 0). In the presence of a magnetic field the non-relativistic charged boson 
gas was studied by Schafroth ||. The relativistic system has recently been extensively 
considered by Daicic et al. ||. However, we do not agree with their conclusions about the 
relativistic Meissner-Ochsenfeld effect. 

According to Ref. [0], introducing an external field is equivalent to introducing an external 
current independent of the dynamics of the system considered. Including the term £ ext = 
Je xt A u in the classical Lagrangian for scalar QED, Euler-Lagrange's equation of motion 
for A v reads 

ff t F ia , = j„+fi* , (1) 

where j u is the induced current. When quantizing the particles but not the gauge field we 
may, in the uniform magnetostatic case, perform the functional integral over the scalar 
field, and write the effective Lagrangian density as 

£eff = £o + £cff C + ^eff + £ext , (2) 

where Cq = —\B 2 is the tree level term, C^ is the one-loop vacuum correction, and C^ 
is the thermal contribution. In terms of the (average microscopic) magnetic induction 
B = V x A, and the external magnetic field H, such that V x H = j ext , we may (upon 
neglecting a surface term) write C cxt = B • H. The mean-field equation then follows from 
a minimization of the effective action (in natural units a = e 2 /47r) 

B = H + M(B) , (3) 

where the average microscopic magnetization is defined by 

w-mD^+tin • (4) 

such that the expectation value of the induced current is (j) = V x M. Considering the 
external field as the acting field (i.e. the field felt by the particles) as in Ref. |I], means 
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that the induced current is neglected in Eq.([L|), and lead to erroneous conclusions about 
the relativistic Meissner effect. In the case of Meissner effect, the magnetization appears 
to be stronger than the external field, as explained below, and thus certainly may not be 
neglected. In the original work by Schafroth || and also in Ref.0, B was used as the 
acting field, consistent with Eqs. (|||4]) . The contribution from the vacuum polarization 
£gj| c to the magnetization is negligible for small magnetic fields, and of no importance 
when considering the topics discussed here, so it will be neglected in what follows. 

In Ref.jT] the thermal part of the effective Lagrangian C^ for the relativistic charged 
boson gas in a homogeneous magnetic field, neglecting all boundary effects, was calculated 
and found to be related to the free energy F, and grand partition function Z(B,T,p,), 
according to C^g = log Z/V/3 = —F/V. Generalizing to d + 1 dimensional spacetime ||, 
the free energy density fd = F d jV d is 

''=-2«^r<i)£f ^(n(EM) + r B (E n m) , (5) 

where f% are the one-particle distributions /g(u;) = (e^™) — l) -1 , and the energy is 
given by 



E n (p) = ^m? + p 2 + (2n + l)eB . (6) 

Following the steps in Ref.0, separating the particle(+) and antiparticle(— ) contribu- 
tions, and introducing dimensionless quantities (J3 = m/3, p = fi/m, fd = fd/m d+1 , pd = 
Pd/m d , B = eB/m 2 , H = eH/m 2 ), we rewrite Eq.@ as 

1 d+i ~ r°° ds _ ( d+i) , /3 2 k 2 . Bs ±kB 

The integral and sum here are absolutely convergent, but if we expand in powers of B, 
the series is only asymptotic, and as p, assumes its critical value at the lowest energy level 
p c = E (0)/m = Vl + B, the coefficients become divergent. However, we may still in a 
simple manner discuss the analytical behaviour of the free energy and the magnetization, 
using the following inequality 

(l + 2Bs)e- Bs >— ^ 7r >(l + Bs)e~ Ss , (8) 

sinn Bs 



which we write as 



US § (l + c Bs)e- Bs , c e [1,2] , (9) 



sinh Bs 



and similarly for the purpose of the magnetization 

8 BS - n 5„ r l 



m^mi-^'^'- c ' e 'f- 21 • (10) 

Introducing the function (x = J3ks) 

^( x )^P[^( x - Xo f + (p cT p)] , Xo = i/(2js c ) , (11) 

and identifying a Jonquiere's function with an exponential argument ip a (z) — SfcLi e~ kz /k a , 
the charge density ep = e(p + + p~) = — e-gL, and the magnetization M d = —Aira-^- are 
written 



t d $ ±(—)^pf^x-^{4,^l 1 ± (x)}+c BPxij^l 1 ± (x)]} , (12) 



l , d+i - r°° dx 
>4nP 

Mf $ -4nacM^)^P 2 f ^-^^^hHx)] • (13) 

The charge density naturally splits into two parts p = p reg + p div , where p reg , the first term 
in Eg. flT2|) , for a vanishing magnetic field is the charge density of noncondensed states, 
and pdiv is the second term in Eq. fll2|) . We shall now investigate pdiv for p — » p c , in 
order to see if the magnetized Bose gas can form a condensate. The leading behaviour of 
■0^3 (7 + ) close to x = Xq is for ^- < 1, or d even || 

d-5 

^^h ± (x)}~T(^)hi^(x-x ) 2 + (p c TP)}\ ' • (14) 

For positive chemical potential (p < can be treated similarly) we thus see that pddn 
remains finite, whereas at p — p c , 4> ±zAll + (x)} ~ (x — x ) d ~ 5 so that p ddiv diverges for 
d < 4. Actually, also p dreg diverges for d < 2. At finite magnetic field, for an arbitrary 
density pd<± there is thus a value of the chemical potential p for any inverse temperature 
P, such that pd = pd,rcg+pd,div, hence the magnetized Bose gas does not condense for d < 4 
(and in the absence of the field the gas does not condense for d < 2). Actually this can 
be seen in a physically more illuminating way. The divergent contribution to the charge 
density should come from the lowest Landau level (n — 0). Separating out that level we 
have according to Eq.([|), after a change of variables of summation and integration 

-+ & f°° J d-3 1 

Pd = - ow-9 w/Qxvd TT / dxx 



2 d - 2 7T d / 2 r(f - 1) Jo exp[p(Vl + x 2 + B - p)] - I 

y / dxx d ~ 3 f!5) 

2<*"%^r(f - 1) fo J v explP(Jl + x2 + (2n + 3)B - p)] - 1 



The first term is easily seen to diverge at pi — p c for d < 4, and the second term is 
exactly of the same form as if the lowest Landau level was included, but with the lowest 
energy (yl + 3B) always larger than the chemical potential (fi < p c ), for finite B. It thus 
follows from the above inequalities that this sum always is finite. Hence even though no 
true Bose-Einstein condensate can form for d < 4, the lowest Landau level can play the 
role of the groundstate and accommodate a large charge density 

Let us now consider the physically most relevant case of d — 3, suppress the dimensional 
index and return to dimensionful quantities. In the case of vanishing magnetic field 
p div = 0, and p rcg is finite as p — > p c (B = 0) = m. At fixed temperature there is thus a 
critical density 

Tm 2 °° 1 
p c (T) = p r egO = m, T, B = 0) = — — Y^ T K 2(kf3m) smh(k(3m) . (16) 

71 fc=i K 

Condensation occurs if this critical density is exceeded. In the limit [i — > /i c the integrals 
in Eqs.(|12],|13|) are dominated by x close to x , and we obtain the leading behaviour 
(fi c = \/m 2 + eB) 



"- s S7f^ BT te ■ (17) 



M S -^ceBL^ 



87rv/2 Hc\] p c - \i 

We now whish to consider the magnetization in the limit of vanishing magnetic field 
B — > 0, at fixed supercritical density p > p c , i.e when a condensate is formed. Then we 
must have that p —>■ p c , and in this limit p rcg — > p c , so that Eqs.(|TT|,|i"BD give that the 
magnetization is approaching a constant. The obtained magnetization law 

M(B-,0) = ~^(p-p e ) , (19) 

2m Cq 

is exactly of the Schafroth form || , who derived it in the non-relativistic case with the 
constant Ci/co = 1. We may actually determine the limiting values of cq and c\ here, 
using the previous method of separating out the lowest Landau level. Taking the limit 
5 — > for all higher Landau levels we find the supercritical density 

eB f°° , 1 

P = W-t I d P 



2n 2 Jo ^ exp[(3(E (p) - p)} - 1 
— J dpp 2 (fB(^p 2 + m 2 ) + f B (^Jp 2 + m 2 



(20) 

fj,=m 



and the magnetization 

M = - e ^r jE - • (2i) 

4W0 E (p)exp[(3(E (p)-p)}-l K ' 

In the n = term we let p approach \Jm 2 + eB in such a way that a prescribed p is 
obtained. We thus find Cq = C\ = 2, in agreement with Shafroth's || non-relativistic 



result. In the ultra-relativistic (T ^> m) and non-relativistic (T <C m) limit, Eq.([16]) 
gives the well-known result 

,J^ 2m ' T>>m (22) 

1 C(|)(^(Tm) 3 / 2 , T«m l } 

Below the critical temperature we can therefore write the magnetization 

e f [1 - (T/T c ) 2 } , T^m 
M(B -»• 0) « p I l v / c; j , 23 

1 ; 2m F [ [l - (T/T c ) 3 / 2 ] , T < m 

where T c is defined by the condition that p c (T = T c ) = p. There is thus a critical field 
H c = —M(B —>■ 0), such that for external fields smaller than H c Eq.(D) has no solution, 
and the field is expelled from the Bose gas. This is the well-known Meissner-Ochsenfeld 
effect. The expulsion of the field is caused by surface currents, and since we do not consider 
surface effects here, they appear as external currents in our formalism. The origin of the 
Meissner-Ochsenfeld effect is here, as well as in the non-relativistic case, that the total 
free energy is minimized if a condensate is formed, which is only possible in vanishing 
magnetic field, and is not connected to the high temperature pair production as claimed 
in Ref. ||. The perfect expulsion of the externally applied field increases the free energy 
per unit volume by || Hl ppl /2 , which causes penetration at supercritical field strengths. 
If we instead keep p = m fixed as B — > 0, Eq.([L8|) gives (p c ~ m + eB /2m) 

M ~ -e—VeBT . (24) 

The square root magnetization law in Eq.(^4j) was obtained in Ref.[|J[], where the corre- 
sponding c\ = 6^71" C[ — 1/2, 1/2] ~ 1.6 was found. Considering the externally applied 
field as the acting field made the authors of Ref. draw erroneous conclusions about the 
Meissner effect from Eg. fl2"4|) . We claim that the magnetization law of Eq. (|I~9D , leading to 



^^Here an expansion of the form [i = m + g{j3, f3 c , B), for g <C eB /m was used. In the limit B — > this 
is only valid for g = 0, and thus this corresponds to keeping the chemical potential (and not the density 
as was intended) fixed for small magnetic fields. 



genuine Meissner-Ochsenfeld effect, may be derived from the expressions in Ref.[|j] if a 
supercritical density is correctly held fixed. 

There has recently been some discussion about condensation of K~ mesons in the core of 
neutron stars and its influence on the equation of state [0]. Since there are also very large 
magnetic fields in some neutron stars it is interesting to ask whether the field can influence 
the condensation. If magnetic flux is trapped in the inner core when the protons become 
superconducting it is believed that flux tubes are formed with field strengths much higher 
than on the surface. The details of the formation and dynamics of such flux tubes are 
not known so we shall only estimate some relevant quantities related to the kaons. We 
can here safely put T = compared to the kaon effective mass m* K ~ 210 MeV and use 
Eq.(^3l) to compute the critical field strength. For the typical value px — 0.1 fm -3 we 
find H c ~ 10 12 T which is far above the maximal observed fields on the surface of neutron 
stars, but comparable with the field strength of the flux tubes in the proton condensate. 
The kaons are essentially non-relativistic so we can use Schafroth's || formula for the 
penetration depth d = [rriK / ^^Pk) 1 — 3 fm. Let us estimate the order of magnitude 
of the field gradient in the flux tube by H c /d. Such a field gradient exerts an enormous 
force even on small magnetic dipoles. For instance, we estimate the force on the electron 
and the neutron, due to their intrinsic magnetic moments p e and p n , to be O(10 3 ) N and 
(9(1) N, respectively. 
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